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Abstract 

We report on two-particle charge-dependent correlations in pp, p-Pb, and Pb-Pb collisions as a func¬ 
tion of the pseudorapidity and azimuthal angle difference, Arj and A(p respectively. These correla¬ 
tions are studied using the balance function that probes the charge creation time and the development 
of collectivity in the produced system. The dependence of the balance function on the event multi¬ 
plicity as well as on the trigger and associated particle transverse momentum (pj) in pp, p-Pb, and 
Pb-Pb collisions at y^snn = 7, 5.02, and 2.76 TeV, respectively, are presented. In the low transverse 
momentum region, for 0.2 < pj < 2.0 GeV/c, the balance function becomes narrower in both Atj 
and A^ directions in all three systems for events with higher multiplicity. The experimental findings 
favor models that either incorporate some collective behavior (e.g. AMPT) or different mechanisms 
that lead to effects that resemble collective behavior (e.g. PYTHIAS with color reconnection). For 
higher values of transverse momenta the balance function becomes even narrower but exhibits no 
multiplicity dependence, indicating that the observed narrowing with increasing multiplicity at low 
Pj is a feature of bulk particle production. 


(c) 2015 CERN for the benefit of the AEICE Collaboration. 

Reproduction of this article or parts of it is allowed as specified in fhe CC-BY-4.0 license. 


See Appendix lAl for the list of collaboration members 





Charge-dependent correlations in pp, p-Pb, and Pb-Pb 


ALICE Collaboration 


1 Introduction 


Angular correlations between two particles have been established as a powerful tool to study the proper¬ 
ties of the system created in high energy collisions of hadrons and nuclei EM. These measurements 
are usually performed in a two dimensional space as a function of At] and Acp. Here At] and A(p are the 


differences in pseudorapidity rj = —In 


tan(0/2) 


(where 6 is the polar angle of a particle relative to the 
beam axis) and in azimuthal angle (p of the two particles. 

In heavy-ion collisions at both the Relativistic Heavy Ion Collider (RHIC) HO] and at the Large 
Hadron Collider (LHC) 1 12 - l3l , these correlations exhibit characteristic structures: (a) a peak at (At] ,A(p) 
= (0,0), usually referred to as the near-side jet peak, resulting from intra-jet correlations as well as cor¬ 
relation due to decay of resonances and quantum statistics correlations, (b) an elongated structure over 
Arj at A(p = 7i originating partially from correlations between particles from back-to-back jets and from 
collective effects such as anisotropic flow, and (c) a similar component at A(p = 0 extending to large val¬ 
ues of Arj, usually called the near-side ridge, whose origin was subject of a theoretical debate ilTHfllll . 
Although initially the near-side ridge was also attributed to jet-medium interactions il7l420ll. it is now 
believed to be associated to the development of collective motion 124143 ill and to initial state density fluc¬ 
tuations, including the initial state effects within the framework of the Color Glass Condensate (CGC) 

iiH-lzl]- 


Similar structures have recently been reported in two-particle correlation analyses in smaller systems. 
In particular, the CMS Collaboration, by studying angular correlations between two particles in Arj and 
A(p, reported the development of an enhancement of correlations on the near-side (i.e. A(p = 0) in high- 
compared to low-multiplicity pp collisions at ^/s = 7 TeV that persists over large values of Arj 0211 . 
In the subsequent data taking periods at the LHC, similar ridge structures were observed on both the 
near- and the aw ay-side in high-multiplicity p-Pb collisions at ^/^nn = 5.02 TeV i33l438n. The origin 
of these effects, appearing in small systems, is still debated theoretically. In particular, it was suggested 
in 091 - 14 in that in high-multiplicity collisions the small system develops collective motion during a short 
hydrodynamic expansion phase. On the other hand, in i42l-l44l1 the authors suggested that the ridge 
structure can be understood within the CGC framework. 

The ALICE Collaboration also reported a particle mass ordering in the extracted V 2 (i.e. the second 
coefficient of the Louder expansion of the azimuthal distribution of particles relative to the symmetry 
plane) values for K^, and p(p) in high-multiplicity p-Pb collisions ll45n . This mass ordering be¬ 
comes evident once the correlations observed in the lowest multiplicity class are subtracted from the 
ones recorded in the highest multiplicity class. The ordering is less pronounced, yet still present, if this 
subtraction procedure is not applied. Similar mass ordering in Pb-Pb collisions idol is usually attributed 
to the interplay between radial and elliptic flow induced by the collective motion of the system. These 
observations in p-Pb collisions were reproduced by models incorporating a hydrodynamic expansion of 
the system 1471. l48n. Recently, it was suggested in 14911 that the signatures of collective effects observed 
in experiments could be partially described by models that couple the hot QCD matter created in these 
small systems, described as an ensemble of non-interacting particles, to a late stage hadronic cascade 
model. More recently, the CMS Collaboration demonstrated that the effects responsible for the observed 
correlations in high-multiplicity p-Pb events are of multiparticle nature iSOll . This strengthens the picture 
of the development of collective effects even in these small systems. 

The c harg e-dependent part of two-particle correlations is traditionally studied with the balance function 
(BP) is in . described in detail in Section |4] Such studies have emerged as a powerful tool to probe the 
properties of the system created in high energy collisions. Particle production is governed by conserva¬ 
tion laws, such as local charge conservation. The latter ensures that each charged particle is balanced 
by an oppositely-charged partner, created at the same location in space and time. The BP reflects the 
distribution of balancing charges in momentum space. It is argued to be a sensitive probe of both the 
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time when charges are created |M, 52] and of the collective motion of the system |^, 53]. In particular, 



In this article, we extend the previous measurements 15411 by reporting results on the balance function in 
pp, p-Pb, and Pb-Pb collisions at .^^nn = 7, 5.02, and 2.76 TeV, respectively. The data were recorded 
with the ALICE detector 155145711 . The results are presented as a function of multiplicity and transverse 
momentum (px) to investigate potential scaling properties and similarities or differences between the 
three systems. The article is organized as follows: Section |2| briefly describes the experimental setup, 
while details about the data sample and the selection criteria are introduced in Section [3] In Section |4| 
the analysis technique and the applied corrections are illustrated. In Section [5l the specifics about the 
estimation of the systematic uncertainties are described. Section 0 discusses the results followed by a 
detailed comparison with models to investigate the influence of different mechanisms (e.g. unrelated 
to hydrodynamic effects) on the balance functions. In the same section, the comparison of the results 
among the three systems is presented. 


2 Experimental setup 

ALICE Q is one of the four major detectors at the LHC. It is designed to efficiently reconstruct and 
identify particles in the high-paiticle density environment of central Pb-Pb collisions |5^, 5^. The ex¬ 
periment consists of a number of central barrel detectors positioned inside a solenoidal magnet providing 
a 0.5 T field parallel to the beam direction, and a set of forward detectors. The central detector systems 
of ALICE provide full azimuthal coverage for track reconstruction within a pseudorapidity window of 
|t] I < 0.9. The experimental setup is also optimized to provide good momentum solution (about 1% at 
Pi <\ GeV/c) and particle identification (PID) over a broad momentum range 


Lor this analysis, charged particles were reconstructed using the Time Projection Chamber (TPC) 16111 
and the Inner Tracking System (ITS) l57ll . The TPC is the main tracking detector of the central barrel 
I 6lll . consisting of 159 pad rows grouped into 18 sectors that cover the full azimuth within |t] | < 0.9. The 
inner and outer radii of the detector are 85 and 247 cm, respectively. The ITS consists of six layers of 
silicon detectors employing three different technologies. The two innermost layers, positioned at r = 3.9 
and 7.6 cm, are Silicon Pixel Detectors (SPD), followed by two layers of Silicon Drift Detectors (SDD) 
at r = 15 cm and 23.9 cm. Linally, the two outermost layers are double-sided Silicon Strip Detectors 
(SSD) at r = 38 cm and 43 cm. 


A set of forward detectors, the VO scintillator arrays |]62]1 . were used in the trigger logic and the mul¬ 
tiplicity determination. The VO consists of two systems, the VOA and the VOC, positioned on both 
sides of the interaction point along the beam. They cover the pseudorapidity ranges 2.8 < T] < 5.1 and 
—3.7 < T] < —1.7 for the VOA and the VOC, respectively. 


Lor more details on the ALICE detector setup and its performance in the LHC run 1, see |]57l. l60ll. 


3 Analysis details 

This analysis is based on data from pp, p-Pb, and Pb-Pb collisions. The data were recorded for pp col¬ 
lisions during the 2010 run at ^/s = 7 TeV, for p-Pb collisions during the 2013 run at .^snn = 5.02 TeV, 
and for Pb-Pb collisions during the 2010 and 2011 runs at .^/^nn = 2.76 TeV. In p-Pb collisions, the 
nucleon-nucleon centre-of-mass system was shifted with respect to the ALICE laboratory system by a 
rapidity of -0.465 in the direction of the proton beam. Lor simplicity, the pseudorapidity in the laboratory 
frame is denoted, throughout this article, with rj for all systems (note that for pp and Pb-Pb collisions 
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the laboratory and the centre-of-mass systems coincide). 


Minimum-bias p-Pb and Pb-Pb events were triggered by the coincidence between signals from the two 
sides of the VO detector. For the pp run, the minimum-bias bigger definition was modified to require 
at least one hit in the SPD or either of the VO detectors. In addition, for Pb-Pb, an online selection 
based on the VO detectors was used to increase the number of events with high multiplicity. An offline 
event selection exploiting the signal arrival time in VOA and VOC, with a 1 ns resolution, was used to 
discriminate background (e.g. beam-gas) from collision events. This led to a reduction of background 
events in the analyzed samples to a negligible fraction (< 0.1%) for all systems BfiOll . All events retained 
in the analysis had a reconsbucted primary vertex position along the beam axis {Zvtx) within 10 cm from 
the nominal interaction point. Finally, events with multipk reconsbucted vertices were rejected, leading 
to a negligible amount of pile-up events for all systems ifiOn . 


After all the selection criteria, approximately 240 x 10^ 100 X 10^ and 35 X 10^ events were analyzed 
for pp, p-Pb, and Pb-Pb, respectively. 


Tracks are reconstructed from a collection of space points (clusters) inside the TPC. The tracking al¬ 
gorithm, based on the Kalman filter, provides the quality of the fit by calculating its value. Each 
space-point is reconsbucted at one of the TPC padrows, where the deposited ionazation energy is also 
measured. The specific ionization energy loss {{AE/Ax)) is estimated by averaging this ionization over 
all clusters associated to the track. The procedure has an uncertainty, which we later refer to as 


To select primary tracks with high efficiency and to minimize the conbihution from background tracks 
(i.e. secondary particles originating either from weak decays or from the interaction of particles with 
the detector material), all selected backs were required to have at least 70 reconsbucted space points 
out of the maximum of 159 possible in the TPC. In addition, the x^ P^r degree of freedom per TPC 
space point of the momentum fit was required to be below 2. To further reduce the contamination from 
background tracks, only tracks with a distance of closest approach (DCA) to the primary vertex in both 
the xy-plane (DCAxy) and the z coordinate (DCA^) below a thi'eshold value (i.e. DCA^y < 2.4 cm and 
DCAz <3.0 cm) were analyzed. These requirements lead to a reconsbuction efficiency of about 80% 
for primary particles and a contamination from secondaries of about 5% at px = 1 GeV/c 16311 in pp 
collisions. The efficiency is similar in p-Pb collisions and it is lower by about 3-5% in central Pb- 
Pb collisions, according to detailed Monte Carlo simulations. In addition, elecbons originating from 
/-conversion and TT^’-Dalitz decays were removed based on the energy loss {AE/Ax) measured by the 
TPC. Tracks for which the measured AE/Ax lied within 30 ^of the Bethe-Bloch paramebization of 
{AE/Ax) for electrons and at least 3a(j£away from the relevant paramebizations for pions, kaons, and 
protons, were removed. 


All particles were reconstructed within |t]| < 0.8. This selection excludes possible biases from the 
backing efficiency that becomes lower for |t]| > 0.8 as compared to |t]| < 0.8. The particles selected in 
this analysis have a bansverse momentum in the range 0.2 < px < 15.0 GeV/c. 


In order to reduce the contribution from track splitting (i.e. incorrect reconsbuction of a signal produced 
by one back as two tracks) and merging (i.e. two nearby tracks being reconsbucted as one back) in the 
active volume of the TPC, a selection based on the closest distance of two backs in the TPC volume was 
applied when forming particle pairs. This was done by excluding pairs with a minimum pseudorapidity 
difference of |At]| < 0.02 and angular distance |A(p*| < 0.02 rad. Here A(p* is the angular distance 
between two backs, accounting also for their curvature due to their charge, according to: 


A(p* = (pi-(p2-ai + a2, (1) 

where (pi and (p 2 are the azimuthal angles of the two tracks at the vertex, and a,- (with / = 1,2) is given 
by 
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CC; — (Ji 


arcsin 


0.0075g^(T)r(cm) 

PT/(GeV/c) 


( 2 ) 


In Eq. |2j qi and q 2 stand for the charge of each track, is the magnetic field in the z direction, r 
corresponds to the radius of the smallest distance of the tracks in the detector used (0.8 < r < 2.5 m with 
a step of Ar = 0.2 cm, for the TPC) and pxi and pt 2 are the transverse momentum values of the two 
particles forming the pair. 

3.1 Multiplicity classes in pp, p-Pb, and Pb-Pb collisions 

The analyzed events were divided into multiplicity classes using the VOA detector. Since this detector 
does not provide any tracking information, the amplitude of the signal from each cell, which is propor¬ 
tional to the number of particles that hit a cell, was used as a proxy for multiplicity lli^l . The choice of 
the VOA as the default multiplicity estimator was driven by the fact that in p-Pb collision^ this detector 
is located in the direction of the Pb-ion and thus is sensitive to its fragmentation li^l . In addition, this 
choice allowed for reducing autocorrelation biases introduced when the multiplicity class was estimated 
in the same rj range as the one used to measure correlations. Eor consistency, the same multiplicity 
estimator was used for the other two systems. Eor the VO detectors, a calibration procedure lOOL 16211 (i.e. 
gain equalization) was performed to account for fluctuations induced by the hardware performance, and 
for the different conditions of the EHC machine for each running period. 

Eor each multiplicity class, the raw transverse momentum spectrum for charged particles with pj > 0.2 
GeV/c reconstructed in | T] | <0.8 was extracted. These raw spectra were corrected for detector a ccep tance 
and efficiency using Monte Carlo simulations with PYTHIA |^], DPMJET |^], and HIJING |£7] event 
generators for pp, p-Pb, and Pb-Pb, respectively. The AEICE detector response for these events was 
determined using a GEANT3 ll68n simulation. In addition to the reconstruction efficiency, a correction 
related to the contamination from secondaries originating from weak decays and from the interaction 
of particles with the material of the detector was applied. This correction was estimated with both the 
aforementioned simulations and also using a data-driven method, based on fitting the DCA distributions 
with templates extracted from Monte Carlo for primary particles and secondaries originating either from 
weak decays or from the interaction of other particles with the detector material, as described in lifil] . The 
resulting corrected charged-paiticle multiplicity was calculated by integrating the corrected transverse 
momentum spectrum over the region with pj > 0.2 GeV/c. 

Table [U presents the multiplicity classes in terms of percentage of the multiplicity distribution, and the 
corresponding number of charged particles with pj> 0.2 GeV/c reconstructed at |t]| < 0.8 for all three 
systems. The resulting values for Ncharged are subject to an overall tracking efficiency uncertainty of 
4% S. 


4 Balance function 


The charge-dependent correlations are studied using the balance function 15111 for pairs of charged par¬ 
ticles with angular differences At] and Acp. Eor each pair, the first (“trigger”) particle has a transverse 
momentum px.trig, while the second (“associated”) charged particle has a transverse momentum px, assoc- 


The associated yield per trigger particle is then calculated for different charge combinations. Eor one 
charge combination (-t,-), it is defined as 

'Note that ALICE also recorded Pb-p collisions but this sample was smaller than the one analysed and reported in this 
article. 
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Multiplicity classes 

PP 

(Ncharged) (coiTected) 
P-Pb 

Pb-Pb 

70-80% 

4.1 ±0.2 

11.2±0.4 

45±2 

60-70% 

5.0±0.2 

16.3 ±0.7 

103 ±4 

50-60% 

6.1±0.3 

18.5±0.7 

204 ±8 

40-50% 

7.4±0.3 

24.1 ±1.0 

364±15 

30-40% 

9.0±0.4 

29.0 ±1.2 

603 ±24 

20-30% 

11.0±0.4 

34.7 ±1.4 

943 ±38 

10-20% 

13.8±0.6 

41.9±1.7 

1419±57 

0-10% 

18.7±0.8 

56.3 ±2.3 

- 

5-10% 

- 

- 

1918±77 

0-5% 

- 

- 

2373 ±95 


Table 1: Corrected mean charged particle multiplicities (for pj> 0.2 GeV/c, and |t]| < 0.8) for event classes 
defined by the percentage of the VOA multiplicity distribution for pp, p-Pb, and Pb-Pb collisions at = 7, 

5.02, and 2.76 TeV, respectively. 


_ 1 d^AassOC,— _ „ Ir 

^(+-) - dAridAq) ~ (+-)//(-r-) 


(3) 


and similarly for the other charge combinations. The signal _) = l/Atrig,+d^Asanie,(+,-)/dAT]dA(p is 
constructed from the number of positive trigger particles Atrig,+ and the particle pair distribution 
d^Asame.(+.-)/dAT]dA<p, formed in Arj-Aq) with positive and negative particles from the same event. Both 
terms are corrected for detector inefficiencies and contamination from secondary particles on a track-by- 
track basis, using the corrections described in Section l3TT] as an inverse weight. is computed after 

summing separately over all events the two components Atrig,+ and d^Asame,(+,-)/d'^^dA(p. 

The background distribution /(+.-) = o:d^A),,jrerf,+-/dAT]dA(p corrects for particle pair-acceptance. It is 
constructed by combining a trigger particle from one event with associated particles from other events. 
This procedure is known as the event mixing technique. These mixed pairs are formed from events 
having the same multiplicity classes and Zvtx within ±2 cm of each other. Each trigger particle is mixed 
with associated particles from at least 5 events. The coefficient a in Eq.[3]is used to normalize the mixed- 
event distribution to unity in the At] region of maximal pair acceptance. Einally, the associated yield per 
trigger particle is computed by calculating the weighted-average of the corresponding yields for several 
intervals of I 4 . This is done to account for the different pair acceptance and efficiency as a function of 

' 4 - 


The balance function is then defined as the difference of the associated yields per trigger particle for 
unlike and like-sign combinations iSlH . according to 


B{Ari,A(p) 




(4) 


The resulting two-dimensional distributions are projected separately onto Arj and A<p and the widths, 
aAj 7 and CJAcp, are calculated as the standard deviation of the distributions. In this analysis, the projection 
in Arj is done on the near- (—71/2 < A(p < TzjT) and on the away-side ( 7 r /2 < Arp < 37 r/ 2 ), separately. 

Three transverse momentum intervals are used in the analysis: the low (0.2 < pt, assoc < Prarig < 2.0 GeV/c), 
intermediate (2.0 < px, assoc < 3.0 < px.trig < 4.0 GeV/c), and high (3.0 < px, assoc < 8-0 < Px,trig < 
15.0 GeV/c) Px regions. Note that the integral of the balance function reported in this article does not 
reach unity but rather 0.5 due to the requirement imposed on the px of the “trigger” and the “associated” 
particles. 
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For 0.2 < pt, assoc < Px.trig < 2.0 GeV/c, the width in At] and A(p is calculated in |At]| <1.6 and — 7r/2 < 
Acp < njl. For higher values of transverse momentum, the balance function distributions are fitted with 
a sum of a Gaussian and a constant. The width is then calculated within 3aGauss> with aoauss extracted 
from the Gaussian of the aforementioned fit. The statistical error of the width is calculated using the 
subsample method 0,|7i. The values of and aA(p are calculated for each subsample (maximum 10 
subsamples were used) and the statistical uncertainty is estimated from the spread of these independent 
results. 


Category 

Systematic uncertainty 

(max. value) 


PP 

p-Pb 

Pb-Pb 

Magnetic field 

- 

- 

1.5% 

FHC periods 

1 .1% 

<0.1% 

1 .0% 

Tracking 

1 .2% 

0 .2% 

1 .2% 

VO equalization 

<0.1% 

- 

- 

Elecfron variafion 

<0.1% 

0 .1% 

0 .2% 

Split/merged pairs variation 

<0.1% 

0 .2% 

0.7% 

Efficiency and contamination correction 

0.4% 

0.4% 

1 .1% 


Table 2: The maximum value of the systematic uncertainties on the width of the balance function over all multi¬ 
plicity classes for each of the sources studied for the pp, p-Pb and Pb-Pb systems. 


5 Systematic uncertainty 

In all figures except Fig. [TJ the data points are plotted with their statistical and systematic uncertainties 
indicated by error bars and open boxes around each point, respectively. The systematic uncertainty 
was obtained by varying the event, track, and pair selection criteria, as will be explained in the following 
paragraphs. The contribution of each source was calculated as the spread of the values of each data point, 
extracted from variations of the selection criteria. If statistically significant, each contribution was added 
in quadrature to obtain the final systematic uncertainty. Following this procedure, the resulting maximum 
values of the systematic uncertainty over all multiplicity classes and systems for the balance function 
projections in At] and A<p were less than 5%. In what follows, we report the maximum systematic 
uncertainties over all multiplicity classes for each system for Oatj and (TAcp- 

The Pb-Pb data samples were analyzed separately for two magnetic field configurafions. The difference 
of 1.5% in fhe resulfs was taken as a systematic uncertainty. For all systems, different FHC periods, 
reflecting different machine conditions and detector configurations (e.g. non-working channels), were 
analyzed separately. The corresponding maximum systematic uncertainties over all multiplicity classes 
was 1.1%. Furthermore, the influence on the results of different tracking strategies was studied by 
repeating the analysis using tracks reconstructed by the combination of signals from the TPC and the 
ITS. The relevant maximum systematic uncertainties from this source were 1.2%, 0.2%, and 1.2% for 
pp, p-Pb, and Pb-Pb, respectively. Finally, the contribution coming from the VO gain equalization in pp 
collisions was investigated by equalizing the signal per VO ring, per channel, and per detector. The study 
did not reveal any systematic differences in the obtained results. 

In addition, several of the track quality criteria defined by fhe tracking algorithm described in Section |3] 
were varied. The uncertainty related to the electron rejection criterium was studied by varying the re¬ 
quirement on the expected Bethe-Bloch parameterization of the momentum dependence of dF/dr for 
electrons from 3a to 5a. This contribution was negligible in the pp system, while it was 0.1% and 0.2% 
for p-Pb and Pb-Pb, respectively. The requirement on the closest distance of two tracks of a pair in the 
TPC was varied from At] = 0.01 to At] = 0.03 and from A(p* = 0.01 rad to A<p* = 0.03 rad. This source 
was found to yield negligible systematic uncertainty for the pp system, while the maximum contribution 
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Fig. 1: The balance function B{Arj,A(p) for charged particles with 0.2 < px.assoc < Pr.trig < 2.0 GeV/c in Pb-Pb, 
p-Pb, and pp collisions at = 2.76, 5.02, and 7 TeV, respectively. From top to bottom the 0 — 5% for Pb-Pb 
and 0 — 10% for p-Pb and pp collisions, 30 — 40%, and the 70 — 80% multiplicity classes are shown. 

for p-Pb and Pb-Pb systems were 0.2% and 0.7%, respectively. The systematic uncertainty of the track- 
by-ttack correction for efficiency and contamination was estimated from Monte Carlo simulations. For 
this, the results of the analysis of a sample at the event generator level (i.e. without invoking either the 
detector geometry or the reconstruction algorithm) were compared with the results of the analysis over 
the output of the full reconstruction chain, using the corrections for detector inefficiencies and accep¬ 
tance discussed in Section [3l This source resulted into a partially correlated uncertainty of around 0.4% 
for the case of pp and p-Pb, and 1.1% for the Pb-Pb system. 

The resulting values for the systematics are summarized in Table [2l for all systems. The table provides 
the maximum value for every source over all multiplicity classes and transverse momentum ranges. 

Finally, different multiplicity estimators were used to study the variations coming from the multiplicity 
class definition. There was no systematic uncertainty assigned for this contribution. The results obtained 
with the two forward detectors (e.g. VOA and VOC) show no significant difference. On the other hand, 
a slightly weaker narrowing of the balance function with increasing multiplicity is observed when the 
central barrel detector is used for both measuring the correlations and the multiplicity class definition, in 
the pp and p-Pb systems. These differences are coming from physics processes (e.g. back-to-back jets), 
whose contribution is reduced if one defines multiplicity classes using a detector located further away 
from mid-rapidity. This also justifies the reason why the VOA detector was chosen as the multiplicity 
estimator in this analysis 
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Fig. 2: The balance function for charged particles with 0.2 < px, assoc < Pr.trig < 2.0 GeV/c as a function of Atj on 
the near-side (upper row) and away-side (middle row) and A(p (lower row) in different multiplicity classes of 
Pb-Pb in panels (a), (d) and (g), p-Pb in panels (b), (e) and (h), and pp collisions in panels (c), (f) and (i) at 
lyiNN = 2.76, 5.02, and 7 TeV, respectively. 


6 Results 

6.1 Balance function in the low transverse momentum region 

Figure [T] presents the balance function for charged particles in At] and A(p for three multiplicity classes 
of Pb-Pb, p-Pb, and pp collisions at = 2.76, 5.02, and 7 TeV, respectively. From top to bottom the 
results for the highest (i.e. 0-5% for Pb-Pb collisions and 0-10% for p-Pb and pp collisions), intermediate 
(i.e. 30-40%), and lowest (i.e. 70-80%) multiplicity classes are shown. The trigger and associated 
particles are selected from the low transverse momentum region 0.2 < px,assoc < PT.trig < 2.0 GeV/c. 
The bulk of the charge-dependent correlation yield is located on the near-side {—7lj2 < Acp < %/2). In 
this region, the balance function becomes narrower with increasing multiplicity for all three collision 
systems. The peak values of the balance function also change with multiplicity, with higher values 
corresponding to collisions with higher multiplicity. On the away-side (7r/2 < Aip < 37r/2), the balance 
function has a larger magnitude for lower multiplicity events. In addition, a depletion in the correlation 
pattern around (At], A<p) = (0,0) starts to emerge in mid-central (e.g. 30-40% multiplicity class) events in 
Pb-Pb collisions and becomes more pronounced in p-Pb and pp collisions with decreasing multiplicity. 
The origin of this structure will be discussed later. 
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Fig. 3: The balance function for charged particles with 0.2 < px, assoc < Pr.trig < 2.0 GeV/c as a function of Arj on 
the near-side (upper row) and away-side (middle row) and as a function of (lower row) for Pb-Pb (panels (a), 
(d) and (g)), p-Pb (panels (b), (e) and (h)) and pp collisions (panels (c), (f) and (i)) compared with results from 
various event generators. Only the highest multiplicity class is shown, i.e. 0-5% for Pb-Pb and 0-10% for p-Pb 
and pp collisions. 


The integral of the balance function over the acceptance is related to measures of charge fluctuations as 
argued in i52l] . and is between 0.25 and 0.35 (i.e. 0.5 and 0.7 in case the pj requirement between the 
“trigger” and the “associated” particles is not imposed) for all systems and multiplicity classes. For each 
system it reveals a mild multiplicity class dependence which, for Pb-Pb, could explain the increase of 
multiplicity fluctuations for central compared to peripheral events reported in iVSil . 


6.1.1 Balance function projections 

Figure |2] presents for Pb-Pb, p-Pb, and pp collisions the projections of the two-dimensional balance 
function in Arj on the near-side (panels (a), (b), (c) ) and away-side (in panels (d), (e), (f)), and A(p in 
panels (g), (h), (i), respectively. The statistical uncertainty, usually smaller than the marker size, is repre¬ 
sented by the error bar while the systematic uncertainty, calculated as the quadratic sum of the correlated 
and the uncorrelated part, by the box around each data point. The balance function as a function of the 
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relative pseudorapidity difference At] on the near-side exhibits a strong multiplicity dependence for all 
collision systems. In particular, the distribution narrows and the peak value becomes larger for high- 
compared to low-multiplicity events. As a function of the relative azimuthal angle Acp on the near-side, 
the balance function exhibits the same qualitative features as for At], i.e. narrower distributions with 
larger magnitude for increasing event multiplicity in all three systems. However, the magnitude of the 
balance function on the away-side exhibits a different trend, with larger values of B(At]) and B(A(p) 
measured for low- compared to high-multiplicity events. 

As aheady discussed in Section [3l in p-Pb collisions, the nucleon-nucleon centre-of-mass system shifts 
by a rapidity of -0.465 with respect to the ALICE laboratory system in the direction of the proton 
beam. The influence of this shift was studied with simulations and, although the balance function is 
not translationally-invariant, the shift does not lead to any significant difference in either the projections 
of the balance function or the extracted widths. 


As indicated previously, starting from mid-central events in Pb-Pb collisions a distinct depletion is ob¬ 
served in the two-dimensional distribution around (At],A<p) = (0,0) that becomes more pronounced in 
events with low multiplicities, and in particular in p-Pb and pp collisions. The fact that the aforemen¬ 
tioned depletion does not seem to be restricted to a very narrow window in either Atj or Acp (the structures 
extend to —0.4 < At] < 0.4 and —7ij6< A(p <7i/6) indicates that the origin is not due to detector effects, 
as was confirmed by independent studies involving modification of cuts controlling track splitting and 
merging. One possible mechanism that could create such a structure is the charge-dependent short-range 
correlations such as Coulomb attraction and repulsion, or quantum statistics correlations 174147611. To test 
this hypothesis, a criterium on the minimum transverse momentum difference Apj between two particles 
of a pair was applied. The value was varied from Apj = 0 GeV/c to Apx = 0.2 GeV/c. The choice 
for the selected values is driven by the fact that the bulk of short-range correlations are expected to have 
Apj <0.1 GeV/c 11771] . The depletion is less pronounced with increasing value of Apj and vanishes for 
Apj = 0.2 GeV/c. The disappearance of the depletion was also achieved by increasing the lower trans¬ 
verse momentum threshold for both the trigger and the associated particle to pj > 0.5 GeV/c. Both these 
observations are inline with the hypothesis that the depletion originates from (mainly) quantum statistics 
correlations and Coulomb effects. The physics conclusion, i.e. narrower distributions with increasing 
event multiplicity, does not change applying one of these criteria. 


6.1.2 Comparison with models 


In Eigs.[3](a),(d) and (g) the balance function in At] on the near- (a) and away-side (d), and in Acp (g) are 
compared with Monte Carlo calculations using the HIJING 16711 and AMPT 178117911 event generators. 
The figures show the 0-5% multiplicity class of Pb-Pb collisions. In AMPT simulations, the string 
melting option was used, with pai'ameters tuned to describe the experimental data on anisotropic flow at 
LHC energies M]. The centrality classes were defined based on the module of the impact parameter. 
It is seen that neither AMPT nor HIJING are able to describe the balance function projections in At] 
on the near-side (see Pig. [3]-a), since they expect not only much broader distributions but they also 
underestimate the magnitude of the balance function. On the other hand, the projection of the balance 
function in At] on the away-side (Pig. |3]-d) indicates that AMPT is in qualitative agreement with the 
data points, contrary to HIJING that predicts a significantly larger magnitude of the balance function. 
Pinally, the A(p projection of the balance function in Pig. [3] (g) shows that while HIJING is still not 
able to describe the data points, AMPT predicts narrower distributions on the near-side but with a much 
smaller magnitude than the one experimentally measured. 

The comparison of the experimental results for the 0-10% multiplicity class in p-Pb collisions with 
model predictions is presented in Fjgs. |3](b),(e) and (h). Por this comparison, results from Monte Carlo 
calculations using the DPMJET and AMPT |78, 79] event generators were used. DPMJET is a 


model based on independent pp collisions, describing hard processes, hadron-hadron interactions, and 
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hadronic interactions involving photons, without any collective effects. This model fails to describe the 
experimental data points in either of the two projections in At], i.e. on the near- and the away-side in 
Fig. [3] (b) and (e), respectively, expecting much broader distributions with smaller (larger) magnitude 
on the near-(away-)side. In addition, for the balance function projection in A(p presented in Fig. |3](h), 
DPMJET predicts broader distributions with a smaller magnitude compared to the measured data points 
on the near-side, but also exhibits a correlation peak on the away-side contrary to what is observed 
experimentally. On the other hand, AMPT, as in the case of the Pb-Pb collisions, seems to describe 
better the balance function projections in both At] and A<p. 


For pp collisions, the experimental results are compared with two variants of calculations using PYTHIAS 
tune 4C 18211 in Figs.[3](c),(f) and (i). This tune contains modified multi-parton interaction (MPI) param¬ 
eters that allow it to describe the multiplicity dependence of (px) ll63n . The default calculation includes 
the color reconnection mechanism, which is switched off in the second configuration. The version of 
PYTHIAS without the inclusion of color reconnection expects a broader balance function near-side pro¬ 
jection in Arj with a smaller magnitude than the one measured as observed in Fig. |3](c). On the other 
hand, the same tune predicts larger magnitude than the one measured for the balance function away-side 
projection in Arj (see Fig.[3](f)). Finally, for the projection in A(p, this tune expects significantly broader 
distributions on the near-side than the measured ones, with an extra correlation peak developing on the 
away-side which is not observed experimentally. On the other hand, the tune of PYTHIAS with the 
inclusion of color reconnection describes the experimental measurement fairly well in both Arj and A(p 
projections. 


As discussed in the previous paragraphs, there are models that exhibit a correlation peak on the away-side 
contrary to what is supported by the data. For this reason, the width of the balance function distribution 
in Arj and A(p will be extracted and compared with models on the near-side only. 


6.1.3 Balance function width 

To quantify the narrowing of the balance function width as a function of multiplicity, the standard devia¬ 
tion a is calculated as described in SectionlH The panels (a), (c), and (e) of Fig.|4]present the evolution 
of Oati on the near-side with multiplicity class, expressed by the multiplicity percentile for Pb-Pb, p-Pb, 
and pp collisions, respectively. Note that the multiplicity decreases from left to right along the horizontal 
axis. The statistical uncertainties of the data points are represented by the error bars and are usually 
smaller than the marker size. For all collision systems, a significant narrowing of the balance function in 
At] with increasing multiplicity is observed. 

The panels (b), (d), and (f) of Figs. |4] show the relative decrease of cJat], expressed by the ratio of 
Cat] for each multiplicity class over the value in the lowest multiplicity class, i.e. 70-80% for all collision 
systems. The narrowing of the balance function with increasing multiplicity is most prominent in Pb- 
Pb collisions where the relative decrease between the largest and lowest multiplicity class is 21.2 ± 
2.4(stat.) zb 2.4(syst.)%. A significant relative decrease is also observed for the other two systems with 
values of 6.7 zb0.2(stat.) zb 0.4(syst.)% and 7.0 zb 0.3(stat.) zb 1.4(syst.)% in p-Pb and pp collisions, 
respectively. Note though that the multiplicities in these thi'ee systems are significantly different (see e.g. 

Table [D 

In Fig. |4](a) the width in At] for Pb-Pb collisions is compared with the results from HHING and AMPT. 
Neither model describes the experimentally observed narrowing of the balance function with increasing 
multiplicity. This is also reflected in Fig.|4](b) where the relative decrease for both models is around 4%. 

Figure |4] (c) shows the comparison of cTati in p-Pb collisions with model calculations. It is seen that 
DPMJET results in broader balance function distributions compared to AMPT. In addition, both models 
expect narrower balance function distributions compared to experimental measurements for low multi¬ 
plicity classes (starting from 60% for DPMJET and 40% for AMPT). However, with increasing multiplic- 


12 





Charge-dependent correlations in pp, p-Pb, and Pb-Pb 


ALICE Collaboration 






p-Pb = 5.02 TeV 
0.75^ • ALICE (c) 

_Systematic uncertainties 

® .DPMJET 

0.65h —AMPT 

0 20 40 60 80 

Muitipiicity ciass (%) 


1.1 


1 - 


° S' 


0.9 


0 . 8 - 


0.2 < p ,p < 2.0 GeV/c 

T,assoc T.tng 

(d) 




0 20 40 60 80 

Muitipiicity ciass (%) 


0.8 


0.75 


0.7 


0.65 


pp \fs = 7 TeV 
• ALICE (e) 

Systematic uncertainties 


0 


0 ..0 




0 


All"" 






A'" 


,X'"' 


PYTHIA 8 

.CR on 

— CRoff 


20 40 60 

Muitipiicity ciass (%) 


80 


1.1 


° S' 


0.9- 


0.8 


0.2 < p < p < 2.0 GeV/c 

T,assoc T,tng 

(f) 





0 20 40 60 80 

Muitipiicity ciass (%) 


Fig. 4: The multiplicity-class dependence of <7^7] in Pb-Pb, p-Pb, and pp collisions at i/s^jn = 2.76, 5.02, and 
7 TeV compared with results from various event generators in panels (a), (c), and (e). Panels (b), (d), and (f) 
show the relative decrease of CTat] calculated with respect to a function of the multiplicity class. The 

transverse momentum values for both the trigger and the associated particles satisfy the condition 0.2 < px.assoc < 
PT.trig < 2.0 GeV/c. 
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Fig. 5 : The multiplicity-class dependence of OAip in Pb-Pb, p-Pb, and pp collisions at i/s^jn = 2.76, 5.02, and 
7 TeV compared with results from various event generators in panels (a), (c), and (e). Panels (b), (d), and (f) 
show the relative decrease of (JA(j) calculated with respect to as a function of the multiplicity class. The 

transverse momentum values for both the trigger and the associated particles satisfy the condition 0.2 < px.assoc < 
PT.trig < 2.0 GeV/c. 
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ity (i.e. below 60% for DPMJET and 30% for AMPT) the balance function distributions are significantly 
narrower in the experiment compared to either of the models. Similar to the Pb-Pb case, neither of the 
models is able to reproduce the significant decrease of the width with increasing multiplicity observed in 
data. This is also reflected in Pig.|4](d), where the relative decrease of the width between the highest and 
lowest multiplicity class for DPMJET and AMPT is marginal and not larger than 2%. 


The experimental results for pp collisions are compared with model predictions in Pig.|4](e). PYTHIAS 
without color reconnection, represented by the solid line, fails to describe the significant narrowing of 
the balance function with increasing multiplicity. The values of Gatj for this calculation are comparable 
within uncertainties to the ones obtained for the lowest multiplicity class in data. On the other hand, 
the inclusion of color reconnection, see the dashed line in Pig. |4] (e), results in a qualitatively similar 
narrowing as the one observed in the measurements. The absolute value of Gat] is lower than the ex¬ 
perimental results for almost all multiplicity classes. Quantum statistics correlations are not included in 
the simulation, which might be the reason for this difference. Pigure |4] (f) that presents the relative de¬ 
crease of Gat] quantifies the previous observations. It is seen that PYTHIAS without color reconnection 
shows a rather weak (i.e. around 2%) narrowing of the balance function with increasing multiplicity. 
This narrowing may result from the increased resonance yield for high- compared to low-multiplicity 
pp events 15411 . The version of PYTHIAS with the inclusion of color reconnection expects a relative 
reduction of around 7%, in quantitative agreement with the measurement. 


Pigure [5] presents the multiplicity dependence of Ga(p in Pb-Pb, p-Pb, and pp collisions in panels (a), 
(c), and (e), respectively. All three systems exhibit a significant multiplicity-dependent narrowing of 
the balance function in Acp. Panels (b), (d), and (f) quantify this nan'owing by presenting the decrease 
of the width in Acp for each multiplicity class relative to the lowest multiplicity class. The data exhibit 
a narrowing of 26.5± 1.0(stat.) ± 1.4(syst.)%, 10.2±0.3(stat.) ±0.2(syst.)%, and lO.S±0.4(stat.) ± 
1.4(syst.)% in Pb-Pb, p-Pb, and pp collisions. 


The multiplicity dependence of the width in Acp in Pb-Pb collisions is compared with expectations from 
HHING and AMPT in Pig.|5](a). HHING fails to describe the experimental measurements while AMPT 
expects a significant decrease of GAcp with increasing multiplicity. The relative decrease in AMPT is 
about 18%, see Pig. [5](b), and can be attributed to a rather strong multiplicity-dependent radial flow in 
the model that acts over the balancing pairs, retaining their initial correlations in Acp. 

The measurements in p-Pb collisions are compared with the results from DPMJET and AMPT in Pig. |5] 
(c). Neither DPMJET, which does not exhibit a significant dependence on the event multiplicity, nor 
AMPT, which exhibits a relative decrease of around 4%, can quantitatively describe the experimental 
findings, as demonstrated in Pig.[5](d). 

Pinally, the values of Ga(p in pp collisions are compared in Pig. |5](e) with the two variants of PYTHIAS 
calculations described before. Similarly to the picture that emerged from the comparison of Gat], the 
variant of PYTHIAS calculation without the inclusion of color reconnection does not describe the strong 
multiplicity dependence reported in pp collision data. However, the calculation with color reconnection 
exhibits a qualitatively similar decrease of Ga(p with increasing multiplicity. The relative decrease for 
this model is around 10%, in quantitative agreement with the experimental results, as indicated in Pig.|5] 

(f). 


The comparison between the data and the corresponding expectations from models like PYTHIA, illus¬ 
trates the potentially significant role of color reconnection on charge-dependent correlations for small 
systems such as pp collisions. The effect of color reconnection in PYTHIAS is strongly connected to 
MPIs, whose number increases with increasing multiplicity. In high-multiplicity pp events, MPIs lead 
to many color strings that will overlap in physical space. Within PYTHIAS approach, these strings are 
given a probability to be reconnected and hence hadronize not independently, but rather in a process 
that resembles collective final-state effects. This process results in a transverse boost of the fragments 
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that leads to the development of final-state correlations between charged particles in a similar way as a 
collective radial boost does. 

6.2 Balance function at high transverse momentum 

In order to study if the narrowing of the balance function is restricted to the bulk particle production at 
low pj or is also connected to hard processes, the balance function was also measured in all collision 
systems for higher values of transverse momentum for both trigger and associated particles. Figure 0 
presents the projections of the two- dimensional balance functions in At] in panels (a), (b), (c), and 
A(p in panels (d), (e), (f) for 2.0 < px, assoc < 3.0 < px.trig < 4.0 GeV/c in Pb-Pb, p-Pb, and pp collisions, 
respectively. The analysis of Pb-Pb and p-Pb collisions was also extended to higher transverse momenta, 
3.0 < Pt, assoc < 8.0 < px.trig < 15.0 GeV/c , shown in Fig.|7] 
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Fig. 6: The balance function for charged particles with 2.0 < px, assoc < 3.0 < px.trig < 4.0 GeV/c as a function 
of Atj (upper row) and (lower row) in different multiplicity classes of Pb-Pb, in panels (a) and (d), p-Pb, in 
panels (b) and (e), and pp collisions, in panels (c) and (f), at = 2.76, 5.02, and 7 TeV, respectively. 

The charge-dependent correlations exhibit little if any multiplicity dependence, in contrast to the findings 
from the lower transverse momentum region. In addition, the distributions in the intermediate and high- 
px range are significantly narrower than the corresponding distributions at lower values of px for each 
multiplicity class. 

The widths of the balance function, Gat] and GAcp for the different transverse momentum regions, are 
presented in Figs. [8] - 19] as a function of the multiplicity class, for Pb-Pb and p-Pb collisions, respec¬ 
tively. The observed narrowing of the balance function with increasing multiplicity is restricted to the 
lower transverse momentum region, i.e. where the bulk of particles are produced. For higher transverse 
momenta, the multiplicity class dependence is significantly reduced, or even vanishes. In addition, the 
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values of cTat] and cJa(p decrease with increasing px for a given multiplicity class. This decrease can be 
attributed to the transition to a region where initial hard-scattering processes and parton fragmentation 
become the dominant particle production mechanism. The emerging hadrons are thus correlated within 
a cone whose angular size decreases with increasing pj. 

For pp collisions, the widths of the balance function (Jat] and Ga(p are compared with results from 
PYTHIA in Fig. [TOl The tune of PYTHIAS without the inclusion of color reconnection is found to 
describe the data at a quantitative level, for both Gatj and GAcp- On the other hand, PYTHIAS with the 
inclusion of color reconnection shows a broadening of the distributions with increasing multiplicity in 
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Fig. 11: The width of the balance function in Atj (a) and in (b) for the three systems analyzed (pp, p-Pb, 
and Pb-Pb), as a function of the charged-particle multiplicity, estimated with the "VOA for |t 7| < 0.8 and px > 
0.2 Ge'V/c. The low-, intermediate-, and high-px intervals correspond to 0.2 < px,a.s.soc < Px.trig < 2.0 Ge'V/c, 
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both At] and Atp, which is not supported by the data. 

6.3 Comparison between the three systems 

A comparison of the widths of the balance function in pp, p-Pb, and Pb-Pb as a function of particle 
multiplicity can provide direct information about differences and similarities between these systems in 
e.g. particle production mechanisms. It is important to note though, that this is performed for different 
center-of-mass energies which could complicate the comparison. 

Figure dU presents the charged-particle multiplicity dependence of the width of the balance function in 
Arj (a) and Acp (b) for all three systems. The results for the low-, intermediate- and high-px intervals 
are shown in the same plot. Multiplicity is defined as the number of charged particles reconstructed in 
|t]| < 0.8 and pj > 0.2 GeV/c, as described in Section |3l It is seen that between the pp and the p- 
Pb systems, and for overlapping multiplicities in the low-px region, the width in both Arj and A(p has 
similar values. This could indicate that the charge-dependent correlations have similar origin in these two 
systems. On the other hand, the comparison of the results between p-Pb and Pb-Pb at the overlapping 
multiplicities indicate differences for both Gatj and (to a smaller extent in) Ga(p- The origin of the charge- 
dependent correlations probed with the balance function in Pb-Pb collisions is believed to be related to 
radial flow and/or to a delayed hadronization scenario. The differences observed in the results of the Pb- 
Pb system compared with the ones in pp and p-Pb collisions at similar multiplicities could be explained 
by a different mechanism that drives the charge-dependent correlations in smaller systems. 

With increasing values of transverse momentum, the balance functions become narrower and exhibit 
no significant multiplicity dependence for all systems, as discussed previously. The origin of these 
correlations at these transverse momentum ranges could be connected to initial hard parton scattering 
and subsequent fragmentation. The agreement of the values of both Gaj] and Ga(p for all multiplicities 
over all three systems clearly indicates that the dynamics responsible for the high-px charge-dependent 
correlations do not change significantly between pp, p-Pb, and Pb-Pb. 

The narrowing of the balance function in both Arj and A(p is a distinct characteristic of the low transverse 
momentum region. Figure [12] visually illustrates the relative decrease of the different systems in this 
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Fig. 12: The multiplicity-class dependence of the width of the balance function in Atj (a) and in (b) for the 
three systems analyzed (pp, p-Pb, and Pb-Pb) relative to the 70-80% multiplicity class. 


region. It is interesting that the relative decrease in this representation is similar between the two small 
systems (around 7% and 10.5% in At] and A(p, respectively) while different for Pb-Pb. The results from 
the analysis of Pb-Pb collisions illustrate a significantly larger relative decrease of 21.2 % in At] (26.5% 
for A(p). A direct comparison of the width at the same multiplicity class can not be done because, for the 
same class, the physics conditions are quite different for pp, p-Pb, and Pb-Pb collisions. However, the 
comparison of the relative decrease and the agreement of the results in both At] and A(p between the two 
small systems could indicate that they share a similar mechanism which is responsible for the decrease 
of the width with increasing multiplicity. 


7 Summary 

This article reports the first measurements of the balance function for chai'ged particles in pp, p-Pb, 
and Pb-Pb collisions at the LHC measured with the ALICE detector. Lor all three systems, the balance 
function in both relative pseudorapidity (At]) and relative azimuthal angle (Acp) was studied for up to 9 
multiplicity classes, and different trigger and associated particle transverse momentum. The widths of the 
balance functions in At] and A(p were found to decrease with increasing multiplicity for all systems only 
in the low-pj region (for pj < 2.0 GeV/c). Lor higher values of pT, the multiplicity-class dependence is 
significantly reduced, if not vanished, and the correlations of balancing partners are stronger with respect 
to the low-pT region. Models incorporating collective effects, such as AMPT, reproduce the narrowing 
of the experimental points qualitatively in Acp, but fail to reproduce the dependence in At]. On the 
other hand, models based on independent pp collisions such as DPMJET and HHING do not show any 
narrowing in p-Pb and Pb-Pb. The comparison of the results in pp collisions with different PYTHIAS 
tunes indicates the importance of MPIs and of the color reconnection mechanism, whose inclusion within 
this model allows for a qualitative description of the experimentally measured narrowing with increasing 
multiplicity at low values of transverse momentum. 
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